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I.  INTRODUCTION 


In  Technical  Report  No.  22  we  presented  a  theory  of  the 
propagation  of  finite  amplitude  ultrasonic  waves  In  hexagonal  crystals. 
In  that  report  we  presented  a  unique  representation  of  the  anisotropy 
of  the  nonlinearity  In  the  basal  plane  of  a  hexagonal  crystal  and 
observed  that  pure  mode  propagation  occurs  along  any  direction  In  the 
basal  plane  as  well  as  along  the  symmetry  axis.  However,  In  presenting 
the  expressions  for  the  nonlinearity  parameters  for  wave  propagation 
along  these  directions  we  overlooked  the  fact  that  an  additional  set 
of  directions  exist  In  which  pure  mode  propagation  can  be  observed. 

These  directions  are  tangential  to  a  cone  whose  apex  angle  Is  bisected 
by  the  symmetry  axis,  and  whose  apex  angle  Is  dependent  on  the  second- 
order  elastic  constants.  The  relative  complexity  of  these  statements, 
the  desire  to  present  a  means  of  defining  the  criteria  for  pure  mode 
propagation,  and  the  desire  to  analyze  the  nonlinear  wave  equation 
for  propagation  In  all  pure  mode  directions  in  crystals  of  hexagonal 
symmetry  as  well  as  trigonal  symmetry,  have  necessitated  our  beginning 
this  report  with  the  fundamental  definitions.  In  most  cases  the  same 
fundamental  definitions  used  In  Technical  Report  No.  22.  He  hope  that 
this  repetition  will  add  to  clarity  and  make  the  present  technical 
report  an  Independent  document.  Although  we  do  repeat  most  of  the 
results  of  Technical  Report  No.  22  In  new  tabular  form,  we  use  a  new 
perturbation  approach  to  the  derivation  of  the  nonlinear  terms,  and  we 
do  not  repeat  the  detailed  equations.  Thus,  even  though  we  consider  the 


present  technical  report  a  complete  document,  we  do  not  consider  It  a 
substitute  for  Technical  Report  No.  22  for  scientists  Interested  In  the 
propagation  of  ultrasonic  waves  of  finite  amplitude  In  crystals  of 
hexagonal  symmetry.  The  two  technical  reports,  then,  provide  Informa¬ 
tion  In  sufficient  detail  for  a  complete  understanding  of  nonlinearity 
in  all  pure  mode  directions  In  crystals  of  hexagonal  symmetry  as  well 
as  nonpiezoelectric  pure  mode  directions  In  crystals  of  trigonal 
symmetry. 

We  begin  the  discussion  with  a  generally  applicable  theory, 
one  capable  of  describing  crystals  of  any  symmetry.  As  a  matter  of 
fact,  at  this  point  one  could  Include  virtually  any  physical  process 
such  as  piezoelectricity  or  any  other  process,  by  Including  the 
appropriate  terms  In  the  expression  for  the  strain  energy.  Our  first 
simplification  Is  to  Include  only  mechanical  nonl Inearl ties  In  our 
strain  energy  expression.  This  specializes  our  consideration  to 
nonpiezoelectric  materials.  The  second  simplification  Is  to  specialize 
to  crystals  of  hexagonal  and  trigonal  symmetry.  This  means  that  our 
equations  are  somewhat  more  complicated  than  those  previously  used 
for  cubic  symmetry,  but  they  still  are  not  so  complicated  as  to  be 
unmanageable.  The  purpose  of  the  present  technical  report  Is  to 
develop  the  theory  appropriate  to  this  situation. 
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II.  GENERAL  THEORY  OF  NONLINEAR  WAVE  PROPAGATION  IN  CRYSTALS 

Consider  a  point  P  In  the  medium  with  coordinates  a^(a,b,c)  In 
the  unstrained  state.  Let  P  move  to  P'  with  coordinates  x.j(x,y,z)  In 
the  deformed  state.  The  components  of  the  displacement  can  then  be 
written  as 

U  *  x  -  a 

V  «  y  -  b  (II-l) 

W  *  2  -  C  . 

In  the  Lagrangian  formulation,  the  strain  Is  described  In  the 
Initial,  or  undeformed,  state  and  the  Initial  coordinates  of  the 
material  particle  a^  are  taken  as  Independent  variables.  The  Lagrangian 
formulation  Is  used  exclusively  In  the  theory  described  In  this  technical 
report. 

The  Lagrangian  strain  parameters  which  are  components  of  the 
finite  strain  tensor  are  given  by  Murnaghan  (1)  as  follows: 


n 


\  (J*J  -  fi) 


(II-2) 


faaa  caggsa  g  gaga 


where  J  and  J*  are  the  Jacobian  matrix  and  Its  transpose;  6  Is  the 
unit  matrix.  They  may  be  expressed,  respectively,  as  follows: 


1  +  Ua  Ub 


J  *  V. 


1  +  u 


1  +  V, 


1  +  V, 


1  +  w_ 


1  +  kl 


(II-3) 


(II-4) 


(II-5) 


where  ua  *  fir*  ub  •  If*  etc- 

The  strain  parameters  can  be  written  In  terms  of  displacement 
gradients  If  we  substitute  (II-3)-( II-5)  Into  (II-2).  Expansion  and 
simplification  give  the  following  expressions  for  n^: 

nn  *  ua  +  1  <ua2  +  va2  +  waZ) 

"12  “  "21  "  ?  'Ub  +  \  *  Vb  +  Vb  +  Vb> 


"13  -  "31  •  r  (Uc  +  ",  +  U.Uc  +  V.Jf  ♦  W 
"22  *  vb  +  J  <Ub2  *  vb2  +  “b2) 


(II-6) 
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"23  ■  "32  *  I  <Vc  4  “b  4  Vc  4  VbVc  4  W 

"33  ■  “c  4  7  <Uc2  4  Vc2  4  Wc2>  ' 

Let  #(n)  be  the  strain  energy  per  unit  of  undeformed  volume. 

The  properties  of  the  crystalline  medium  enter  Into  the  theory  through 
the  strain  energy  density  4>(n)  which  can  be  expanded  and  expressed  as 
a  sum  of  terms  of  different  degrees  In  the  elements  of  n  as  follows: 

♦  *  *0  +  *1  4  *2  +  *3  4  4  (H-7) 

In  terms  of  the  elastic  constants,  Eq.  (I 1-7)  can  be  written  as 

♦  *  ♦()  4  klC1jT11j  4  k2C1jkin1jnki 

+  k3Cijk£mnn1jT,k£Tlmn  +  •**  »  (n-8) 

where  the  C's  are  the  elastic  constants  and  kn  (n  >  1,2,3  ...)  Is  a 
constant  factor  depending  on  the  definition  of  elastic  constants.  One 
does  not  lose  generality  by  setting  the  first  two  terms  In  (II-8)  equal 
to  zero  since  4q  Is  the  energy  of  the  undeformed  medium  and  corres¬ 
ponds  to  displacement  without  deformation.  Thus,  the  expression  for 
the  strain  energy  density  becomes: 

♦  *  *2  4  * 3  4 

“  k2C1jkin1jnki  4  k3C1jktmnn1jT,kinmn  *  (IN9) 

As  Is  well  known,  are  the  second-order  elastic  (SOE) 

constants  and  are  the  third-order  elastic  (TOE)  constants. 


fS 
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In  accordance  with  Brugger's  thermodynamic  definition  of  elastic 
constants  (2),  the  factor  kn  is  1/nl.  Therefore,  based  on  that 
definition  we  have  kg  *  ^  and  k^  *  Using  a  contracted-subscript 
notation  for  the  elastic  constants  as  well  as  using  Brugger's 
definition  of  elastic  constants,  the  general  expression  for  <t>  can  be 
rewritten  as: 

♦  ■  1  CiA  A*  +  i  CAyyn1  At^mn  + 

higher  order  terms  ,  (11-10) 

where  a  -*■  1j,  p  ->  kJ.  and  v  -*•  mn.  For  example,  we  write  1  for  the  pair 
of  indices  11,  2  for  22,  3  for  33,  4  for  23  and  32,  5  for  31  and  13,  6 
for  12  and  21.  Summation  over  repeated  indices  is  understood. 

The  SOE  constants  form  a  fourth  rank  tensor  containing  81 
components,  of  which  21  are  Independent  for  the  most  unsyumetrlcal 
triclinic  crystal  and  the  TOE  constants  C>yv  form  a  sixth  rank  tensor 
with  729  components,  of  which  56  are  independent  for  the  trl clinic 
crystal.  The  number  of  elastic  constants  decreases  considerably  for 
crystals  of  higher  symmetry.  Eq.  (11-10)  for  the  strain  energy  density 
takes  the  appropriate  form  for  crystals  of  different  symmetry.  The 
number  of  elastic  constants  In  different  crystal  classes  has  been 
worked  out  by  various  authors  and  has  been  presented  In  a  tabular  form 
by  Heannon  (3). 

The  equations  of  motion  for  an  elastic  medium,  in  fact,  are  a 
restatement  of  Newton's  second  law.  For  convenience  one  can  Introduce 
the  stress  tensor  T,  which  Is  not  symmetric,  as 
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T  -  <> 


(n-n) 


where 


3  <t>  „ 
3n 


3$/  3Hj  -j 

34>/3n-j  £ 

34>/ Sn-j  ^ 

34>/3^21 

34>/3n22 

3$/  3t^2  3 

3  4>/  3n£2 

34»/3n33 

(11-12) 


Thus,  the  equations  of  motion  In  Lagrangian  coordinates  are 
written  as  (4): 


3Tij 

aa.  =  p0Ui  * 

J 


(II-13) 


These  equations  of  motion  take  particular  forms  along  the  a,  b  and  c 
axes  of  the  crystal : 


ii  =  hi  +  !!li  +  9T13 

p0  da  3b  3c 
•• 

"0V  *  -5T- +  IT  +  lir  • 

ST^i  ST^p 

*0W  ■  TT  +  TT  *  1? 


(11-14) 


According  to  Eq.  (11-11),  the  stress  matrix  T  can  be  written  as: 


11 

T12 

T13 

jn 

J12 

J13 

Bf 

*»ii 

3<f 

an,  j 

3<p 

3n13 

T  „ 

t 

m 

j 

J 

ii 

21 

'22 

23 

J21 

u22 

u23 

*"21 

3n22 

3n23 

31 

T32 

T33 

J31 

J32 

J33 

3n31 

3  4> 

3n32 

3f 

8n33 

(11-15) 


From  (11-15)  we  can  write  down  the  components  as: 
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T11  '  J 


T12s  J 


t13  =  J 


'21 


'22 


'23 


'31 


T32  =  J 


T33'  J 


11 

O 

Snl  1 

+ 

J12 

an2, 

+ 

J13 

ay 

n31 

34 

+ 

J 

34 

+ 

.i 

34 

11 

Sn12 

d12 

3n22 

°13 

3n32 

34 

+ 

34 

+ 

J.« 

34 

11 

Sn13 

J12 

3n23 

u13 

3n33 

34* 

+ 

j 

34 

+ 

3<J> 

21 

8nn 

J22 

an2. 

°23 

3n31 

L_ 

94 

+ 

.i 

34 

+ 

.i 

34 

21 

3”12 

d22 

3n22 

J23 

3n32 

L, 

9  4 

+ 

.1 

34 

+ 

.i 

34 

21 

3nn 

u22 

3n23 

J23 

8n33 

L, 

d  $ 

+ 

J  „ 

34 

+ 

34 

31 

an,, 

°32 

an2. 

u33 

an3, 

L, 

a* 

+ 

3  4 

+ 

,1 

34 

31 

an, 2 

°32 

3n22 

u33 

9n32 

L, 

34 

+ 

J  « 

34 

+ 

j 

34 

31 

S"13 

°32 

3n23 

d33 

9n33 

(11-16) 


The  components  of  the  Jacobin  matrix  in  terms  of  displacement 
gradients  can  be  expressed  from  ( I 1-3) : 


Jjl  *  1  +  Ua  *  1  +  3U/3a 
J-12  r  Ub  *  3U/3& 

J13  =  Uc  *  3U/9c 
J21  *  Vfl  *  3V/3a 

J22  *  1  +  Vb  *  1  +  3V/3b 
°23  *  vc  *  aV/3c 


(II-17) 
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*  B  7  **11^11  +  n??  +  nl?  +  ^91  )  +  ^12^nnn22  ~  7  n12^  ”  5"  n2^ 


+  C13*nlln33  +  n22n33^  +  7  C33n332  +  Wn132  +  n312  +  n232  +  n322^ 
+  1  Cni(nll  +  n22)3  +  1  C113n33(r)n  +  n22)Z  +  7  C133n332(nn  +  n22^ 

+  C144(n11  +  n22)(n232  +  n322  +  n]32  +  n312) 


+  ?  cififi^nn  +  )C (^n  “  +  2(r,i22  +  n212^ 


(III-l) 


+  6  C266^nll  “  n22^  “  6^n122  +  n212^nll  "  n22^ 

1  3  2  2  2  2 
+  6  C333n33  +  C344n33*n23  +  n32  +  n13  +  n31  ^ 

+  7  C366T)33^nll  ‘  n22^2  +  2^n122  +  n212^ 

2  2  2  2 

+  2C456^n22  “  nll^n23  +  n32  '  n13  "  n31  }  +  4^n12n23n31  +  n21n32n13^ 


In  this  expression  Cjggt  C2gg,  C3gg  and  C4gg  are  combinations  of  TOE 
constants  given  by 


C1 66  =  T  (‘2Cm  ■  C112  +  3C222  ^ 
C266  "  1  (2C111  "  C112  *  C222^ 


C366  *  7  ^C113  "  C123^ 


7  (_C144  +  C155*  * 


(UI-2) 


Elnspruch  and  Manning  (8)  derived  an  expression  by  using  Birch's 
definition  of  TOE  constants.  That  expression  can  easily  be  converted 


Basal  Plane 


Figure  III-l.  Hexagonal  close  packed  structure. 
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These  equations  are  to  be  examined  so  one  can  define  the  directions 
and  evaluate  nonlinearity  parameters  for  finite  amplitude  ultrasonic 
waves  propagating  in  hexagonal  and  trigonal  crystals. 


III.  PURE  MODE  PROPAGATION  IN  HEXAGONAL  CRYSTALS 


There  are  two  classes  of  crystals  with  hexagonal  symmetry,  one 
with  the  Hermann-Mauguin  symbol  6,  <T,  6/m  with  5  second-order  and  12 
third-order  elastic  constants  and  the  other  with  the  symbol  622,  6  mm, 
6m2,  6/irmm  with  5  second-order  and  10  third-order  elastic  constants. 

Since  most  of  the  crystals  with  hexagonal  symmetry  belong  to  the  second 
class.  We  will  confine  our  attention  to  It.  The  hexagonal  close  packed 
structure  Is  shown  In  Figure  III-l.  A  drawing  of  the  coordinate  frame 
and  the  basal  plane  are  also  given  In  the  figure.  The  five  SOE  constants 
are  C^,  C^,  C^,  C33  and  C44  and  the  ten  TOE  constants  are  C^,  C^3, 

C116*  C1 33*  C144*  C226*  C333*  C334*  C366  and  C456  By  us1n9 

Brugger's  definition  of  third-order  elastic  constants,  the  expression 

for  the  strain  energy  density  for  hexagonal  crystals  Is  given  by  (7): 


where  ^  and  *3  represent  the  separate  parts  of  the  strain  energy 
density  containing  only  SOE  constants  and  only  TOE  constants, 
respectively. 

Now  as  soon  as  the  expression  of  the  strain  energy  density  in 
terms  of  strain  components  Is  acquired  for  specific  crystals  under 
Investigation,  we  will  be  able  to  get  the  wave  equations  for  propagating 
along  any  a"-direct1on  In  the  crystals  by  using  the  above  equations. 

These  expressions  of  the  strain  energy  are  written  out  in  more 
explicit  terms  for  crystals  belonging  to  hexagonal  and  trigonal 
symmetries  in  the  next  sections. 

Combining  (11-37),  (11-38)  and  (11-39),  we  may  have  the  following 
stress  components  in  which  the  strain  terms  involving  SOE  constants  are 
separated  from  those  involving  TOE  constants: 


«*4>2  3^2 

+ 


*♦2  +  2  av"  +  a*3 

anjji  an; ^  aa "  ^21 


(11-40) 


a*2  a*2  ar  a*3 

an^  ^^11  "  ^^31 


Therefore,  we  may  write  a  set  of  equations  of  motion  for  plane  sound 
waves  propagating  along  the  a"-axis  direction  in  crystals: 
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Moreover,  the  components  of  the  Jacobian  matrix  Is  terms  of 
displacement  gradients  can  be  reduced  to 


1  +  3U"/3a" 
aV'/aa" 
aw"/ aa" 


(II-37) 


J12  *  J13  *  J23  *  J32 


Thus,  the  stress  components  Involved  in  Equations  (11-19)  also  can  be 
reduced  to: 


Tn  *  Jn 

T21  *  °21  8^  +  3^7 
T31  '  J31  8^  +  8^  ’ 


(11-38) 


Obviously,  under  the  above  rotated  coordinate  system,  the  strain 
energy  density  $  can  be  expressed  as  a  function  not  only  of  the  double 
primed  strain  components  n]|j  but  also  the  rotation  angles  $  and  e;  i.e., 
It  can  be  written  as: 


♦(nij)  *  ♦(njj.t.e) 


*2^n1 *  ♦g(n^j»$»6)  > 


(11-39) 
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nll  *  ^.|Cos2teos2e  -  (n'^'3  +  n3-j)cos2*cosesine 

‘  ^n12  +  n21^cos<,sin*cose 
o 

n12  *  nVlcos4>s1n^cos  6  *  W3  +  n3-j  )cos4iSln4>coses1n6 
+  (n^2cos24.  -  n21s1n24>)cose 


n13 


2  2 

T^'-|Cos<l>co$esine  +  (n^cos  e  -  n3-|Sln  e)cos$ 
-  n^sln^slne 


2 

n21  “  cos  ♦sin  ♦cos'  e  -  (n}'3  +  n31)s1n^cos^s1necose 

2  2 

-  nV2s1n  Kose  +  nj^cos  ♦cose 

2  2  2 
n22  *  niisin  <|,cos  6  "  ^n?3  +  ♦slnecose 

+  (n}'2  +  n21)s1n^coSKOse 

2  2 

n23  *  -|sin*coses“,ne  +  (nf3cos  e  -  n3-|Sln  e)s1n$ 

♦  n21cos^sine 

2  2 

n31  *  *  (n][3sin  e  -  n31cos  s)cos$ 

-  n^sln^slne 

2  2 

n32  *  T»iis1ni*'s',necose  *  (n}'3s1n  e  -  n31cos  e)s1n$ 

♦  n^cos^stne 

2 

n33  “  6  +  ^ni3  +  ^31  Jsinecose  . 


(11-36) 
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Thus,  Eq.  (11-32)  can  be  simplified  as: 


8  cos e(r^-j cose  -  ursine)  -  n^sinecose 


n12  “  ni2cose 


n13  *  (T»i'is'ine  +  n^cose)  -  Ti“1s1n2e 


n21  *  n2lcose 


"22  -  0 


(11-35) 


n23  *  ^i5"*06 

2 

*  (n^cose  -  n}'3s1ne)s1ne  +  n31cos  e 
n32  * 

n33  *  (n^slne  +  n}’3cose)sine  +  n^cosesine 

Putting  (11-35)  Into  (11-31),  we  have  the  following  expressions 
for  the  strain  components  In  the  original  frame  expressed  in  the  terms 
of  the  strain  components  In  the  twice-rotated  or  double-primed  coordinate 
system: 


t  V 
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n23  *  sin$  nj3  +  cos<j>  n23 
n3j  *  n3j  cos$  *  n32s1n$ 

n32  *  n31  s’n*  +  t132cos^’ 

n33  =  n33 

and 

*  cosetn^cose  -  n^slne)  -  s1ne(n31cose  -  n33s1n0) 
n|2  *  n^cose  -  n32s1ne 

ni3  *  <”n  sine  +  n^'3cose)cose  -  (n^slne  +  n33cose)sine 
n21  *  n2lcose  *  n23s1ne 

nj2  *  n22  (11-32) 

n23  *  n2is,*ne  +  n23cos 6 

n31  *  (I)iicose  *  n][3s1ne)s1ne  +  (n31cos0  -  n33s1ne)cose 
n32  *  r‘i2s1ne  +  n32cose 

n33  *  (n^slne  +  r^'3cose)s1ne  +  (n31s1ne  +  n33cose)cose  . 

If  only  the  propagation  of  the  plane  wave  along  the  a"-ax1s  Is 
accounted  for,  the  displacement  components  will  be  a  function  of  only 
the  coordinate  an  (see  Eq.  (11-27)).  In  this  case,  we  will  have: 

n22  *  n33  *  n23  *  n32  *  0  ,  (11-33) 

and  the  other  nonvanishing  displacement  gradients  will  be  simplified 
as  follows: 


and 


(n)  «  (R^)(n,)(R1) 


(11-28) 


(n* )  -  (R|)(n")(R2) 

to  arrive  at  the  following  expressions: 

cost  -slnt  0  nj-j  nj2  ^3  cos*  sin*  0 

n  *  sint  cost  0  n21  n22  n23  -slnt  cost  0  (11-29) 

0  0  1  n31  n32  n33  0  0  1 

and 


cose 

0 

-sine 

nll 

n12 

-H 

n13 

cose 

0 

sine 

0 

1 

0 

n21 

.H 

n22 

_  II 

n23 

0 

1 

0 

sine 

0 

cose 

_  M 

"31 

.H 

n32 

n33 

-sine 

0 

cose 

(11-30) 

Multiplication  and  simplification  lead  to  the  following 
expressions  for  the  strain  components  expressed  In  terms  of  the  strain 
components  In  the  twice-rotated  system. 

n-j-j  *  cost (n{j cost  -  nj2s1nt)  -  s1nt(n21cost  -  n22s1nt) 
n-|2  *  cost(n]-|Slnt  ♦  nj2cost)  -  s1nt(n2-jSlnt  ♦  n22cost) 
n13  *  njjcost  -  n23s1nt 

n21  ■  s1nt(njjCost  *  nj2s1nt)  ♦  cost(n21cost  -  n22s1nt) 


t\99  •  s1nt(n^s1nt  ♦  n^cost)  ♦  cost(n4,s1nt  ♦  n«cost) 


(II-31) 
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Then  we  may  rotate  the  above  rotated  basal  plane  (a'b1)  by  an 
angle  e  about  the  b’  axis  (Figure  11-1)  so  that  the  new  coordinate 
frame  formation  can  be  expressed  as 


(11-24) 


a" 

a* 

b"  -  (R2) 

b' 

c" 

c* 

where  the  rotation  matrix  (R2)  Is  given  by 


cose  o  sine 

0  1  0 

-sine  0  cose  , 


(11-25) 


and  Its  corresponding  transpose  (R|)  Is 


(R|)  *  0 


cose  0  -sine 

0  1  0 

sine  0  cose  . 


(11-26) 


Let  us  define  the  displacement  components  propagating  along  the 
a"-ax1s  In  the  new  coordinate  system  as 


U"  «  U"(a\t) 

V"  -  V"(a",t) 
r  »  W"(a"ft)  . 


(H-27) 


Using  the  transformation  matrix  as  given  by  (11-22) •  (11-23),  (11-25) 
and  (11-26)  one  can  transform  the  strain  components  by  (1) 


n 


The  procedure  for  accomplishing  this  Is  developed  In  the  next  section 
of  this  technical  report. 

We  have  shown  the  strain  energy  density  t  In  the  original  frame 
can  be  expressed  as  a  function  of  strain  components, 

♦  =  ♦(nij)  .  (11-20) 

Now,  suppose  we  want  to  consider  the  propagation  of  the  plane 
finite  amplitude  wave  along  some  arbitrary  direction  Instead  of  the 
a-axis  under  original  consideration.  We  can  rotate  the  basal  plane  In 
the  original  frame  by  an  angle  t  about  the  c-axis  (Figure  II-l)  first 
so  that  the  rotated  coordinate  frame  formation  can  be  written  as 


a* 

a 

b' 

*  (R-j ) 

b 

c' 

c 

where  the  rotation  matrix  (R^)  Is  given  by 


(Rj ) 


cos*  sin*  0 

•sin*  cost  0 

0  0  1 


and  Its  corresponding  transpose  (R|)  Is 


cos  -slnt 
slnt  cost 
0  0 


0 

0 

1 


(11-21) 


(11-22) 


(11-23) 
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that  keeping  higher  order  terms  Is  not  justified.  The  terms  kept  in 
the  theory  to  be  presented  are  sufficient  to  account  for  most  of  the 
observed  nonlinear  properties. 

III.  PURE  MODE  DIRECTIONS 

In  an  anisotropic  medium  there  are  only  certain  directions 
along  which  elastic  waves  can  propagate  as  the  pure  longitudinal  modes 
which  are  of  primary  interest  In  the  application  of  the  harmonic 
generation  technique.  Associated  with  any  propagation  direction  there 
are  possibly  three  independent  waves*  the  displacements  of  which  form 
a  mutually  orthogonal  set  (coupling  terms).  In  general*  none  of  the 
three  displacement  vectors  coincides  with  the  vector  which  is  normal 
to  the  wavefront;  1.e.»  In  general,  the  waves  are  neither  pure 
longitudinal  nor  pure  transverse.  The  specific  directions  In  which 
pure  mode  longitudinal  waves  propagate  have  been  determined  qualita¬ 
tively  by  Borgnis  (5)  and  Brugger  (6)  for  crystals  with  cubic*  hexagonal 
and  trigonal  symmetry*  etc.  In  order  to  derive  the  nonlinear  wave 
equation  of  a  pure  longitudinal  mode  In  specific  directions  and  then  to 
determine  the  third-order  elastic  constants  from  the  harmonic  generation 
technique*  It  is  necessary  to  know  quantitatively  these  pure  mode 
directions.  In  addition,  some  unique  pure  mode  directions  even  depend 
upon  the  SOE  constants  for  the  crystals  under  consideration. 

In  order  to  determine  the  pure  longitudinal  mode  directions, 
one  must  develop  a  way  to  rotate  the  coordinates  and  observe  the 
directions  In  which  the  coupling  terms  In  the  wave  equations  vanish. 


J33  “  1  *  Wc  =  1  +  3W/3C 

We  will  consider  the  case  of  plane  finite  amplitude  waves 
propagating  along  the  axes  of  the  medium  under  consideration.  For 
plane  waves  propagating  along  the  a-axis  the  displacements  become: 

U  *  U(a,t) 

V  *  V(a,t)  (11-18) 

W  *  W(a,t)  . 

As  the  a-axis  can  be  rotated  into  alignment  with  any  axis  of  the  medium, 
the  equations  of  motion  in  general  can  be  written 
•• 

pQU  «  aT^/aa 

P0V  *  aT21/aa  (11-19) 

P0W  =  aT31/aa  , 

for  propagation  along  the  a-direction. 

These  equations  are  written  out  In  more  explicit  terns  and  solved 
for  crystals  belonging  to  hexagonal  and  trigonal  symmetry  in  the  follow¬ 
ing  sections.  The  expressions  of  elastic  strain  energy  density  appro¬ 
priate  to  different  crystal  synmetrles  are  obtained  with  only  terms  up 
to  third-order  since  to  date,  experimental  uncertainty  is  great  enough 


Into  Eq.  (III-l), which  was  Brugger's  definition.  The  relation  between 
Birch  and  Brugger  TOE  constants  can  be  obtained  from  the  relation: 


(Birch)  (Brugger) 

CAUV  *  "W6  (I1I-3) 

where  M  is  the  possible  number  of  ways  In  which  CAyv  can  be  expressed 
in  tensor  notation. 

Now  we  substitute  (11-36)  for  (I1I-1)  to  give  the  expression 
for  the  strain  energy  density  In  the  double  primed  coordinate  system, 

♦  (n^)  *  ^(nj'j,*,0)  +  ♦jlnjj.t.e)  (III-4) 

where 

*2*1  cn^niiC0s4e  +  W2  +  n21  ^cos^® 

-  2n^  (n}'3  +  n^Jsinecos30  +  (n^3  +  n31)2cos20s1n2o] 

-  ^  ci2^ni2  +  n2i^cos2e  +  C13Cnj|2cos20s1n20 

♦  nj'i(nj3  +  nj3)cos0s1n0(cos20  -  s1n2o)  -  (nj3  +  n31)2cos20s1n20] 

+  ^  c33^]is1n4e  +  2nil  (ni3  +  n31)s1n30cos0 

+  (n3  3  +  n31  )2sin20cos2e]  +  C^[2nY2cos2es1n2e 

+  (nj2  +  n!>2)s1n20  +  20^(0^  +  n^jcoseslnefcos2©  -  s1n2o) 

♦  (n.;3  +  ^31  )2(cos4o  sin4©)  -  An^n^cos^sln^]  (III-5] 

and 


,*o  .* 
*  -  . 


*3  *  l  Cni[n^c°s6e  -  3n^(nJ3  +  n31  )cos5es1ne  +  SnfjCnfj  +  n31)Zcos4es1ne 
‘  tn13  +  T>31)3cos3esin3e]  +  ^  Cm[n|3s1n20cos40  +  nj2(nj3  +  n^) 

•  s1necos3e(cos2e -2sin2e)  +  n^(n^3  +  n31)2cos2esin2e(s1n2e  -2cos2e) 

4  (n^3  +  n31)3cos3esin3e]  +  C133[n^s1n4ecos2e  +  n^f  (n^3  +  nj,) 

•  cjsesin3e(2cos2e  -  sin2e)  +  n^(nj3  +  n31)2  cos2es1n2e(cos2e  -2s1n2e) 

-  (n^3  +  n31)3cos3es1n3e3  +  C144[n^3(2cos4esin2e) 

4  nlVn12  +  i2Z)cos2es1n2e  +  2n;2(nJ3  +  n31)cos3esine(cos2e  -  2sin2e) 

‘  (n13  +  T,3i>(T‘ii  +  nj2)cos0s1n30  +  4  nJ2)cos6e 

4  3ni'i(n13  +  r>32)s1n4ecos2e  -  ZnJ^nJ2  4  n32)cos4es1n2e 
4  4nJ1nJ3n31cos2es1r;20{s1n2e  -  2cos2e)  -  (nj3  +  n31) 

•  +  n32)coses1ne{cos4e  4  s1n4e)  +  4nJ3n31(nJ3  +  n31)cos3esin3e] 

4  7  cl66[nllcos6e  ’  3nll(n13  *  n31)cos5esine  4  +  r^)2 

•  cos4esfn2e  -  (nj3  +  n31  )3cos3flsfn3e  4  2nj-|(nj|  4  T>22)cos4e 

-  2(nJ3  4  ^31 ) (nii  +  Ti22)cos3es1ne]  4  C2gg[nJ3cos60cos(24>) 

•  (cos2(2*)  -  3s1n2(2^)) 43nJ2(nJ3  +  n31)cos50sin0cos(2*) 

•  (3s1n3(2$)  -  cos2(2$))  4  +  n31)2cos4o  s1n20cos(2*) 

•  (cos2(2*)  -  3s1n2(2$)  4  3njf (r»J2  4  n51)cos50s1n(2^)(s1n2(2^) 
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-  3cos2(2$))  +  (nl2  +  nJJ^cos^cosU^sIn2^*) 

"  nlVn12  +  r>22)cos4ecos{2<|>)(cos%  +  s1n%)  +  6n^(n\'2  ♦  n21) 

•  ^n13  +  ^31  )c°s4es1nes1n(2<*>) (3cos2(2<».) -s1n2(2*)) 

+  6nJ1ri;2Tj51cos4esin2(2*)cos(2^)  +  (nj3  +  n51)3cos3es1n3ecos(2^) 

•  (cos2(24>)  +  3s1  n (2^> } )  -  (nj2  +  )3cos3es1n3(2$) 

+  3(n^'2  ♦  n2i )(ni3  +  ^3*1  )2cos3es1n2es1n(24>) (s1n2(24>)  -  3cos2(2<^) ) 

-  9(n]|2  +  n2i^n13  +  )cos3es1nes1n2(2<|>)cos (2i*») 

+  ^n13  +  n3l^n12  +  ^l^05305100005^*^005^  +  s1n%) 

+  (nj2  +  n21)(n^'2  +  n22)cos3es1n(2*)(cos%  +  sin4*) 

-  6nf2  n21  ^n13  +  n31  )s,,r,9cos30cos(2*)s1n2(2*) 

•  ®Tt32n21  ^n12  +  n21  )COS^0S^3(2^)]  +  ^  ^333tr'iis^n^®  +  3rjj2 

•  (nj3  ♦  n31)s1n5ecose  +  3*1^  (nf3  +  n31)2  sin4ecos2e 

+  tn13  +  n51)3sin3ecos3e]  +  C344[2n;3cos2es1n4e  ♦  n^tnj2  +  n22) 

•  s1n4e  +  2n^'2(n]|3  +  n31)cos0s1n3e(2cos2e  -  s1n2e) 

♦  nll^n13  +  *»3i)sin2e(cos4e  ♦  s1n4e)  -  4nJ1n;3n31cos2es1n4e 

+  (nj2  +  n22)(nj3  ♦  T»31)s1n3ecose  ♦  2nJi(n}3  +  n31)2sin20cos2o 

•  (cos2e  -  $1n20)  +  (n{3  ♦  r^MnJ2  +  n32)s1necose(cos40  ♦  s1n40) 


♦  n^(nj3  +  n31)cos3es1ne(cos2e  -  2s1n2e)  +  +  n31)2 

*  cos2esin2e(s1n2e  -  2cos2e)  +  (n^3  +  n31)3cos3esin3e 
+  2nJ1(nJ2  +  nJ;2)cos2es1n2e  +  2(nj|  +  r^Mn^  +  )cos3es1ne] 

+  2C456[2n^,3cos4es1n2e  +  2n^'2(n^'3  +  n31  )cos3esine(cos2e  -  2s1n2e) 

-  2n|1  (n|3  +  n31)2cos2es1n2e(cos2e  -  s1n2e)  -  nJjUj2  +  n22) 

•  cos2esin2e(2s1n2(2<(i)  +  cos2(24>))  +  Zi^-jU^'2  -  n32) 

•  coses1n2ecos(24i)s1n(24i)  +  2n^2(n32  -  nj2)coses1n2ecos(24>)s1n(2$) 
+  8nJ1nJ2Ti21cos2esin2e(cos4^  ♦  s1n%  -  s1n24>cos2«») 

+  2n\,1nJ3n31(cos6e  ♦  cos2es1n4e)s1n2(2$)  -  Zn^tn^2  +  nj2) 

•  cos4esin2es1n2(2$)  +  n^Cnj2  +  n32)cos2e(cos4e  +  s1n4e)cos2(2$) 

-  +  n31  )2cos2esin2e(cos2e  -  s1n2e)cos2(24>) 

-  2n^ri2in]f3cos4es1n8cos{2t)s1n(2^)  +  2nJ1n21n31cos2es1n3e 

•  cos(24>)sin(2^)  -  Znj-in^^sInecosScosU^sInU*) 

♦  2n;1n;2n1M3cos20s1n3ecos(2^)sin(2«)  +  2t,;1(t,J2  +  n^)2 

•  cos2es1n2es1n2(24)  -  *ij'|nj3n3jCos4es1n2ecos2(2$) 

♦  (n;2  ♦  +  n22)s1n2ecosecos(2+)s1n(2t) 


-  2(nl|2n21  +  n12T$Vcoses1n  esin(2^)c°s (24»)  +  2(nJ3  +  n31) 

•  (n3|  +  n32)cos3es1n30s1n2(2$)  -  2nJ3n31(nJ3  +  n31) 

•  coses1n0(s1n40  +  cos4e)s1n2(2$)  -  2nj2(nj2  +  n21) 

•  cos3es1n2es1n(2$)cos(2$)  -  (n}'2  ♦  n21  )<nli  *  ^i^056 

•  (cos4e  +  sin4e)cos(24>)sin(24>)  +  4(nj2  +  n21  >^1 3n31 

•  cos30sin2ecos{2<ti)sin(2(|>)  -  (n}'2  +  n21  )n32cos30sin20 

•  s1n(24>)cos(2$)  +  2n^'3n31(n^'2  +  n21)cos0s1ne(cos40  +  s1n40) 

•  sin(2<(»)cos(24>)  +  n21nl  3^n12  +  n2i  )cos3esi  nesi  n2  (2<*) ) 

+  (nf 3  ♦  n3i )cos0s1n30s1n2(2$)  ♦  4nJ2n21(nJ3  +  n31) 

•  cos0sin0(cos20  -  s1n2e)(cos%  ♦  sin%)  -  2(nf2  +  n21 )n21n31 

•  sin30cos0s1n2(2$)  ♦  2n^2n31 (n|2  +  n21 )s1necos3es1n2(2*) 


-  2nj2TiJ3(nJ2  +  n21)s1n30cos0sin2(2«{i) 

-  2(nJ3  ♦  n^HnJ*  +  nJ2)cos3es1nes1n2(2*)  . 


( I I 1-6 


In  order  to  obtain  a  set  of  equations  of  motion  for  the  plane 
waves  propagating  along  the  a"-ax1s  direction,  it  Is  necessary  only  to 
calculate  three  of  the  strain  derivatives  of  the  elastic  energy,  l.e., 
8+/dnfi»  d*/an2j  and  d^/dn3-|t  according  to  (11-38).  Those  strain 
derivatives  can  be  calculated  and  simplified  by  taking  Into  account 

n12  *  n21  and  ni3  and  n3-j.  They  are: 


.  •*'  .  A  'MS1  \  *.*  '  '  •  ’  *  *  «  *  •  * 
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■  _2|_ .  !ig_  ♦  !4- 

j  3nn 


tin-7) 


where 


M  A  A  M 

(C^  cos  0  ♦  2C13  cos^e  s1n*0  +  C33  sin  e 

♦  C44  sin2(2e))n{1  *  (-C^  cos2e  s1n(2e) 

♦  C13  cos(2e)s1n(2e)  +  C33  s1n20  s1n(2e) 


+  C44  s1n(40))n;3 


(II 1-8) 


■gjjir-  *  j  C^^Cn^cos^e  -  An^n^cos^slne  +  4n^'|cos4esine] 
+  j  C113[3n^,2s1n28co$S  +  ^'i^sinocos3© 


•  (cos20  -  2s1n20)  +  4n^'3cos20s1n20(s1n20  -  2cos20] 


4  ?  C133[3n^2s1nScos20+  4n][1n]|3cos0s1n30(2cos20  -  sin2e) 
♦  4n^2cos20s1n20(cos20  -  2s1n20)]  +  C^44[6n][2cos^8s1n20 


♦  2nj2Cos2es1n28  ♦  4n^(nj2  ♦  )cos30sin0(cos20  -  2s1n20) 


♦  2n^2(cos60  ♦  5cos2es1n28cos(20)) 


I 


* m  '*•  A  »V^ mj *  ^ 


m 
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♦  J  ci66^3n?lcos6e  "  6niVni3  +  n31)cos5es1ne 
+  3(nJj  +  n31)2cos4es1n2e  +  2(nj|  +  nJ2)cos4e] 

+  J  C266[nJ2cos6ecos(24')(cosZ(24')  -  3s1n2(24>)) 

+  12n^,1n^3cos5es1necos(24i){3s1n2(24.)  -  cos2 (2^) ) 

+  4n'1'2cos4es1n2ecos(24>)(cos2(2$)  - 3s1n2 (24») ) 

+  12nJ1n^2cos5es1n(24.)(s1n2(24.)  -  3cos2(2$)) 

+  42n^'2cos4ecos(2*)s1n2(24)  +  24nJ2nJ  3cos4es1nes1n(24>) 

•  (3cos2{2$)  -  s1n2(24)]  +  J  C333[n^sin68 

+  4n^n^'3s1n5®cos0  ♦  4nf2s1n4ecos2e] 

2  2  4  3 

♦  C34^[6n^cos  esln  e  +  8n^r>^3coses1n  e 

♦(2cos28  -  s1n2e)  ♦  2n3|s1n^e  +  2nf2(5s1n2ecos4e  +  s1n6e 
-  6cos2es1n4e)]  +  ^  C3gg[3n^2s1n28cos4e 

♦  4^1nJ3s1necos3e(c°s2e  -  2s1n2e)  ♦  4nJ2cos2es1n2e 


vs- 


v'v 


30 


*(s1n2e  -  2cos2e)]  +  2C^gg[6n^cos^esin2e 
+  ®nl 3cos3es‘*ne (cos^0  -  2s1n2e)  +  2n}'2(-4cos4es1n2e 
♦  4cos2es1n*e  +  cos®e  +  5cos2es1n*ecos(2<f>) 


+  cos2es1n*esin2(2<|>)  -  2cos4es1n2esin2(2«ji) 

-  6cos4es1n28cos2(2$) )  +  2n^2COS2esin2e 


•  (sin2{2$)  -  cos2(2$)  +  4cos%  +  4s1n%) 


2  2  2 

+  4n}2n{3COs(2^)s1n(2^)cos  es1ne(s1n  e  -  cos  e)]  , 
^  **z  +  3*3 

8n2i  ^2] 


where 


a*2 

^21 


2  2 

(C^cos  e  -  c^cos  e  +  2C 


44 


sin2e 


(III-9) 

(III-10) 


(HI-11) 


and 


3*3 


2  2  3 

C144t2nll  njjicos  esln  e  -  es*ne^ 

+  F  C266^3nllcos5esin^24,^s<n2^2*^  "  3cos2(2*)) 

+  12nJ1r>][3COS4esines1n(2#)(3cos2(2^)  -  s1n2(24>)) 


•  *  »  *  «  “  •  •  ■  ’  •  ' 

*  ■% *  » 


,*  *.»  '  *  *  •  *  •  "  •  *,*  *  *  *  «  *Ji  *  »  *  »  *,»  '  »  '  ^  *  *  *  »  *  *  *  •  '  »  ^Jj 


38nJinj2c°s4ec°s(2*)s1n2(2$)  -  2n^nj'2COS4ecos(2e) 
•  (cos%  +  s1n^4>)  +  12nj2cos30sin20sin(24>) 


•  (s1n*(2$)  -  3cost{24>)  -  30n^|cos3©sin3(24>) 

+  6n^2COS3es1n{24»)(cos%  +  s1n^4>)  -  76r>i2ni3s^necos3ecos 

•  (24>)sin2(24>)  +  4nY2nY3cos3esfnecos(24»)] 

+  C144[2Ti51nJ2sin4e  +  4n^2n^'3s1n3ecose] 

+  \  C366^nlln12sin2ecos2e  +  ®n12n13s*necos^ 

o  o  9  9  A 

+  2C45g[2n^n^2(cos4es1n^e  +  4cos  esln  ecos % 


♦  4cos20sin2s$in%)  +  2nj’1ny3cos2©sin2ecos(2*)sin(2<j) 
-  3n^2^,3cos3es1nes1n2(2<(i)  -  4n52nJ3coses1n3esin2(24>) 
+  8n^2n^3(cosS  +  s1n%)cosesfne(cos2e  -  s1n2e)] 


3r)31  3n31  3n31 


(III-12) 


(III-13) 


where 


34«  «  22 

9  [-Cj-jcos  eslne  +  C13coses1ne(cos  e  -  sin  e) 
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3  9  9 

♦  C33s1n  ecose  ♦  2C44coses1ne(cos  0  -  sin  0)]^ 

♦  UC^cos^sin2©  -  4C13cos2es1n2e  +  2C33cos2es1n2e 

4  4  22 

♦  4C44(cos e  ♦  sin  0  -  sin  ecos  ®)]n|3 


(II 1-14 ) 


?  cm^"niicos5es1ne  +  ^ni  i  ni  3cos^es1ne 

-  4n^cos3es1n3e]  +  C1^[n^cos3esine(cos2e  -  2s1n2e) 

+  4r)^n^3cos2es1n2e(s1n2e  -  2cos2e)  +  12n^2cos3es1n3e] 

+  1  Cl33^nllcosesin3e^cos2e  ’  s<,nZe) 

+  4n^n^'3cos2es1n2e(cos2e -2s1n2e)  -  12n^,2cos3es1n3e] 

+  C1^[2n^'2cos3es1ne(cos26  -  2s1n2e)  +  2n^ny3cos6e 

+  10n^n^3s1n4ecos2e  -  12n|'1nY3s1n2ecos4e  -  2n^2Coses1n3e 

+  2nf2cos5fls1ne  +  2n},2cosesfn*e  +  72nf3cos3es1n3e] 

1  O  C  i  p 

+  £  C16g[-3n^cos  eslne  +  ^n^n^cos  esln  e 

-  12n^cos3esin3e  -  4n|2cos3es1ne] 

+  ^  C266^3nllcos50s1n6cos^2*^3s1n2^2*^  "  cos^ (24>) ) 

♦  12n^n][3cos^es1n2ecos(2^)(cos2(24>)  -  3sin2(24>)) 

+  12nj|')nj2cos*es*nesM2$)(3cos2(2$)  -  s1n2(2^) ) 

+  12n^2cos3es1n3ecos(2$)(cos2(2$)  +  3s1n2(24>) ) 

-  42n^2cos3es,,nes1n2(24)cos(2$)  -  2nj2cos3es*ne 

•  cos(2$)(cos%  +  s1n%)  +  24n^2T»^'3cos3es1n2es1n(2^) 

•  ($1n2(2*)  -  3cos2(2t))] C33[3nj^sin5e  cose 
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+  12n^nj3sln*ecos2e  +  nj^sln^cos2©] 

+  C344[2nJ2coses1n3e(2cos2e  -  sin2©) 

+  n^'1n^3Sln2e(10cos4e  +  2s1n40  -  12cos2es1n2e) 

+  2r^'2sin3ecose  +  4n^3COses1ne[cos4e  ♦  sin4© 

-  3cos2es1n2e)  +  y  C366^  t fcos3es"*ne (cos^©  "  2s1n2e) 

+  4n^,1n^3cos2es1n2e(s1n20  -  2cos2e) 

+  12nJ3Cos3es1n3e  +  4nj2s1necos3e3 
+  2C456[2n^,2cos3es1ne(cos2e  -  2s1n2e) 

+  n^,1^3(-8cos2es1n2ecos(2e)  +  cos6e  ♦  12cos2es1n4ocos2(2<i>) 
-16cos4es1n2ecos2(24i))  +  n^3Coses1nes1n2(24>)(12cos2es1n2e 

44  o 

-  6s1n  e  -  6cos  e)  +  n^cos (24>)s1n(24>)cos  esl ne 

2  2  3  2 

•  (sin  e  -  cos  e)  +  ^n^n^cos  esln  0cos(24»)s1n(2$) 

2  2  3  2  2  2 

-  2r^2Sln  (2$)s1n  ecos0  +  4n^2coses,lne(cos  0  -  sin  0) 

•  (cos%  +  s1n%)]  (III-15) 

If  we  put  these  expressions  for  the  strain  derivatives  of  elastic 
energy  Into  (11-41),  we  would  get  a  set  of  wave  equations  for  propagating 
along  the  a^-axls  direction.  It  Is  predicted  that  the  set  of  wave 
equations  should  be  nonlinear  and  a  perturbation  approach  Is  necessary 
In  order  to  solve  It. 


Let  us  suppose  the  solutions  of  these  nonlinear  wave  equations 
may  be  expressed  In  the  following  forms: 


U"  -  UJ  ♦  UJ  ♦  ... 

V"  *  V"  +  VJ  +  ...  (III-16) 

W"  «  WJ  +  WJ  +  ... 

where  Uj\  Vj,  Wj'  and  U2»  V2,  Wjj  represent  the  first  and  second 
approximation  solutions  of  the  displacements,  respectively.  First, 
we  look  at  the  first  approximation  equations.  Putting  (III-16)  Into 
(11-41)  and  using  (III-7)-(III-15)  and  (11-34),  we  can  obtain  a  set 
of  first  approximation  wave  equations  for  propagating  along  a"-ax1s 
direction  as  follows: 


p0U?  *  (C11cos4e  ♦  ^  C13s1n2(2e)  ♦  C33s1n40 

32U" 

+  C44s1n2(2e))  — ^  (-C11cos20s1n(2e) 

3a" 

+  C13cos(2e)s1n(2e)  +  C33s1n2es1n(2e)  ♦  2C44s1n(2e)cos(2e)) — 


*o»i  ■  I  <cu 


cos2e  -  C12cos2e  +  2C44s1n2e)  — -j 


p0*ij  *  \  (-C11cos2es1n(2e)  +  C13s1n(2a)cos(20) 

32Um 

♦  C,,s1n20s1n(20)  ♦  2CA«s1n(20)cos(2o))  — i 

"  44  3a"z 

4  ?  ^7  Cnsm2(20)  -  C13s1n2(20)  +  £  C33s1n2(2e)  +  4C44(cos4e 

♦  s1n4e  -  i  s1n2(2e)3 

4 


(III-17) 


which  can  be  rewritten  as: 


32Ui’  32Wi' 

P°U1'“^+62^J 


32v; 

C0V1  -  "1  ^ 


32WV  32UV 

p0^1  *  Y2  ?+  ®2  F 
u  L  3a"^  6  3a 


(III-18) 


Here 

a  *  C11cos4e  +  C13s1n2(2e)  +  C33s1n4e  +  C44s1n2(2e) 

Y1  *  7  (cnC0sZe  *  C12cos2e  +  2C44s1n2e) 

Y2  *  1  C^s1n2(2e)  -  C^3s1n2(2e)  +  ^  C33s1n2(2e) 

+  2C44(1  -  j  sin2 (2e ) ) 

e2  *  ’  ?  C11cos2esin(2e)  +  ^  C13cos(2e)sin(2e) 

+  5-  C33s1n2esin(2e)  +  C44s1n(2e)cos(2e) 


(III-19) 

In  the  above  equations  the  terms  containing  the  coupling 
coefficient  e2  ere  called  the  coupling  terms. 

Obviously,  because  the  coupling  terms  are  existing  In  the  set 
of  wave  equations.  It  becomes  Impossible  to  get  the  solution  of  the 
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pure  longitudinal  vibration  even  when  only  the  pure  longitudinal 
vibration  mode  is  excited  at  the  surface  of  the  transducer. 

In  order  to  seek  the  direction  along  which  the  pure  longitudinal 
mode  can  propagate,  we  seek  the  directions  along  which  the  coupling 
coefficient  e2  vanishes  (or  is  minimum).  It  is  necessary  to  note  that 
since  the  elastic  constants  C^.  are  contained  in  the  expression  for 
the  coupling  coefficient,  the  direction  we  want  to  seek  can  be  deter¬ 
mined  only  for  a  specific  crystal  rather  than  in  general  terms.  But 
the  coupling  coefficient  $2  onty  a  function  of  the  angle  e,  so  it 
means  that  under  the  same  angle  6  the  polar  diagram  of  the  coupling 
coefficient  e2  versus  the  angle  *  should  be  a  circle. 

From  (II 1-19)  we  may  see  Inmedlately  that  $2  vanishes  when 
e  =  m/2  (n  *  0,  1,  2,  ...).  This  reveals  that  the  pure  longitudinal 
mode  can  propagate  along  any  radial  direction  in  the  original  basal 
plane,  as  well  as  in  any  direction  perpendicular  to  the  original  basal 
plane.  Moreover,  the  above  condition  without  coupling  coefficient  $2 
should  work  without  exception  for  any  hexagonal  crystals,  which  Is  In 
coincidence  with  that  In  Technical  Report  No.  22  by  J.  Philip  and 
M.  A.  Breazeale  (7). 

It  Is  Important  to  point  out  that,  besides  the  above  directions, 
another  direction  with  vanishing  e2  also  exists  In  hexagonal  crystals. 

It  Is  the  direction  along  a  cone  whose  apex  angle  -  ep)  is  centered 
about  the  c-axis  (in  original  coordinate  system).  However,  as  mentioned 
above,  the  apex  angle  Is  different  for  the  different  hexagonal  crystals. 
Figure  II 1-2  shows  the  directions  of  pure  modes. 


Figure  IV-1. 


(b)  Set  of  equivalent 

directions  in  a 
trigonot  system 


The  coordinate  system  for  crystals  of  trigonal  symmetry. 
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+  ^n23  +  n32^n122  +  n212^ 


1  2 

+  F  C133n33  ^11  +  n22^  +  C134n33{^nll  *  n22^n23  +  *32* 

+  ^n31  +  n13^n12  +  n21  ^  +  C144  {^n232  +  n322^nll 

2  2  2  2  2  2 
+  n22 ( n31  +  n13  ^  +  C155{n22^n23  +  n32  ^  +  nll^n31  +  n13 


,  1  p  3,1/*  3  ,  /*  //  2  ,  2\ 

6  C222n22  F  C333n33  C344n33Un23  n32  > 

+  (n312  +  n132))  +  ^44^^23  +  n32^  "  ^n23  +  n32^n312  +  n132^} 

1  2 

+  ?  (Cni  +  Cn2  -  C222)n22 

+  F  ^ ~C  114  ‘  2C124^n222^23  +  n32*  +  1  ^"2C111  "  C112  +  3C222^nll 

(n122  +  n212)  +  l(2Cm  '  C1 12  "  C222> 

x  n22(n122  +  n21Z)  +  ^CH3  ”  Cl23^n33^n122  +  n212)  +  F^C114 

+  3^124^11  ^n31  +  ni3^n-j2  +  n21) 

+  7(C114  -  ci24)n22(n31  +  n13)(n12  +  n21}  +  F*'C144  +  C155) 


(n23  +  n32)(n31  +  n13)(n12  +  n21}  + 


(IV-1) 


In  this  expression  ^113’  ^114*  ^124*  ^144  ^222 

are  combinations  of  TOE  constants  given  by 


i *  •  *  -**  »  “  •'*  »  ’  .  •  /•  ,  ■  »"**'•_* •  ,*» V»  , "'»■  „*•  ,*•  .*•  ,*• 
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The  unit  cell  of  a  crystal  with  the  trigonal  symmetry  Is  shown 
In  Figure  IV-1.  A  drawing  of  the  coordinate  frame  Is  also  given  In 
the  figure.  The  six  SOE  constants  are  ,  ^12*  ^13*  ^14*  ^33*  ^44 
and  the  fourteen  TOE  constants  are  ^ i  ^112*  ^113*  ^114*  ^123*  ^124* 

C1 33*  C1 34 *  C144*  C155*  C222’  C333’  C344*  C444* 

Some  of  the  trigonal  crystals  (for  instance,  quartz  and  lithium 

niobate)  exhibit  piezoelectricity  that  affects  their  nonlinear 
behavior,  and  proper  attention  ultimately  must  be  paid  to  it  in  this 
sort  of  crystals  (10),  especially  for  some  strongly  coupled  piezo¬ 
electric  trigonal  crystals  like  LINbOj.  In  this  report,  however,  we 
will  concentrate  on  the  mechanical  properties  and  neglect  the  effect 
of  piezoelectricity.  Piezoelectricity  was  treated  to  some  extent  in 
Technical  Report  No.  23. 

For  the  nonplezoelectrlc  trigonal  crystal,  Kaga  (11)  has  given 
the  elastic  strain  energy  density  with  Brugger  definition  of  third- 
order  elastic  constants.  After  making  some  corrections  (3),  the 
energy  density  can  be  expressed  as. 


♦  *  1  Cll(nll2  +  n22 +  C12nlln22  +  C13*n22n33  +  n33nllJ 
+  C14{(nn  -  n22)(n23  +  n32)  +  (n31  +  n]3)(n12  +  n21)> 

+  \  C33n332  +  C44^n232  +  n322  +  n312  +  n13  * 

♦  \  (Cn  -  Cl2)(n122  +  n212) 

+  F  C111T,113  *  1  C112nll2n22  +  ?  C113(nll2n33  +  n222n33) 

♦  \  C114nll2(n23  +  n32)  4  C123nlln22n33  +  C124{nlln22(n23  +  n32} 


A*. 


*  *  *  '  *  *  *  "  »  ’  »  *  *  ‘  •  '  •  '  »*»*»*  f1*  »  *  k  '  •  »  ,  ^  ^  »  ,**  jJ 
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Figure  1 1 1-5.  Magnitude  of  the  K3  parameter  plotted  as  a  function  of 
angle  In  the  basal  plane. 


.'O  .*  *.  *.  •.  A  ,\ 
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negative  of  the  ratio  of  the  nonlinear  term  to  the  linear  term  In  the 
nonlinear  wave  equation  of  pure  longitudinal  mode  [Eq.  (II 1-34)] .  It 
Is  given  by 

e  *  -  £  .  ( i  n-36 ) 


In  terms  of  A-j  and  A2.  It  can  be  rewritten  as. 


6  = 


(II 1-37 ) 


Therefore,  by  measuring  A-j  and  A2  one  can  determine  B  which  can 
be  used  to  evaluate  the  6's  which  are  combinations  of  SOE  and  TOE  con¬ 
stants,  or  the  K^'s  which  are  combinations  of  TOE  constants  only. 

The  parameters  K 2  and  K3  for  the  specific  directions  of  pure 
longitudinal  mode  also  are  written  In  tabular  form  In  Table  III-l,  In 
which  under  K3  we  Indicate  those  specific  angles  4  which  lead  to  magni¬ 
tudes  of  the  quantity  cos(24)[cos2(24)  -3  sin2  24]  of  1,  0,  or  -1. 

These  would  be  directions  along  which  measurements  will  yield  the  most 
useful  Information.  This  behavior  of  K3  In  the  basal  plane  was  given 
In  graphical  form  In  Figure  2  of  Technical  Report  No.  22,  which  we 
reproduce  here  for  convenience  (Figure  II 1-5) . 


IV.  PURE  MODE  PROPAGATION  IN  TRIGONAL  CRYSTALS 


There  are  also  two  classes  of  crystals  with  trigonal  symmetry, 
one  with  the  Hermann-Mauguln  symbol  3?  and  the  other  with  the  symbol 
32,  3m,  3m.  The  first  class  has  seven  SOE  constants  and  twenty  TOE 
constants,  and  the  second  has  six  SOE  constants  and  fourteen  TOE  con¬ 
stants.  In  this  technical  report  we  limit  our  Investigation  to  the 
latter. 
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a  «  Kg 

and  ( I 11-32 ) 

6  «  OKg  +  K3)  . 

Thus, 

6  2  4  4  2 

K3  *  o  *  C^cos  e  +  3C^13sin4ecos  e  +  3C133sin  ecos  e 

+  12C^4cos4esin2e  +  3C166cos6e  +  C26gCos6ecos(2«j)) 

•  (cos2(2d>)  -  3s1n2(2$))  +  C333s1n^e  +  12C344cos2es1n*e 

+  3C36gSin2ecos*e  +  24C45gCos*es1n2e  .  (111-33) 


Now,  we  may  simplify  Eq.  (11-41)  and,  as  a  matter  of  fact, 
write  the  nonlinear  wave  equation  of  pure  longitudinal  mode  as 


a  32UM  +  6  3U"  #  32U" 
3a^  3a^ 


(III -34 ) 


The  solution  to  the  nonlinear  wave  equation  can  be  expressed  as 


U"  *  U"  +  UJ 

»  AjSlnfk^a"  -  u>t)  -  A2cos2(kLa"  -  wt)  .  ( II 1-35) 

6  2 

Here  A2  «  [j^](kLA.|)  a"  Is  the  amplitude  of  the  generated  second 
harmonic  wave.  The  solution  (III-35)  Is  In  complete  analogy  to  cubic 
crystals;  however,  the  expressions  for  a  and  6,  or  alternatively  Kg 
and  K3  must  be  examined  for  hexagonal  synmetry. 

As  in  the  case  of  cubic  crystals  (7),  we  can  define  the 
ultrasonic  nonlinearity  parameter  for  hexagonal  crystals  as  the 


Now,  equating  the  coefficients  on  both  sides  of  the  equation 


we  get 


-  u,2P()Ba"  +  k20Ba"  +  kLoC  '  -  J  6  k2A2 
and 

-  u>2P0Ca"  -  kLaB  ♦  K^aCa"  =  0  . 

2 

But  a  *  pgC^  and  u>  *  C^k^,  so  Eqs.  (II 1-26)  become 

-  cfkfp0Ba"  +  kfp0cfto"  ♦  k^cjx  ■  -  |  kft 


and 


CLkLp0Ca"  ‘  kLp0CLB  +  kLp0CLCa"  *  0  • 
Dividing  Eq.  (III-27)  by  ^PqC2,  we  get 

fiCk^)2 


P  r  2 


0CL 


Also,  dividing  Eq.  (I 11-28)  by  kpQCL2,  we  get 


B  =  0  . 


( I 11-26 ) 


(III-27) 


(III— 28) 


(II 1-29) 


( I 11—30) 


So,  the  solution  to  the  wave  equation  of  second  approximation 
(II 1-22)  becomes 


«(k.A,)2 

*  -  C — LJ7-]a"cos2(k1aH  -  «t)  (III-31) 

8p0  CL"  L 


If  we  use  the  same  symbols  as  Breazeale  and  Ford  (9)  used  before  for 
cubic  crystals,  we  can  put 


-  jv  ■r:  *\ r.  ^  w.  jv 


T7  ".l-pi'vi  -  t  -  r 


T-T-T 
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a2UJ;  auv  a2uy 

pou2  *  °  — 7  +  6  aT"  — 7 
u  c  3a "c  38  aa"£ 


(III-21 ) 


or 


a2u" 

p0^2  "  a  ~J.  *  "  ?  fikLAls1n2*kLa"  '  wt) 

3d 


(II 1-22) 


Here 

6  *  3a  +  o  (II 1-23 ) 

and 

6  2  4  4  2 

a  *  C^^os  e  +  2C113$1n  ecos  e  +  3C133s1n  ecos  e 

+  12C144cos4esin2e  +  3C166cos6e  +  C2ggcos60cos(2(|») 

•((cos2(2*)  -  3sin2(24>))  +  C333s1n6e  +  12C344cos20s1n40 

♦  3C3g6s1n2ecos4e  +  24C43gCOS4es1n2e  .  ( I 11-24) 


Taking  into  account  the  boundary  condition  for  U£,  that  Is, 

U!J  *  0  where  a"  *  0,  we  use  a  trial  solution  in  the  following  form: 

UJ  «  BaMsin2(kLa"  -  u>t)  +  Ca,,cos2(kLaM  -  at)  (III-25) 


where  B  and  C  are  the  coefficients  to  be  determined.  Putting  the 
solution  Into  (III -22) •  we  get 

4(-w2pQBaH  ♦  k2aBa "  +  kLaC)sin2(kLa"  -  wt) 

«•  4(-«2p0CaM  -  kL«B  ♦  k[aCa")cos2(kLa"  -  u>t) 

2  / 

—  s1n2(kLa"  -  <ut)  . 


7 
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The  coefficient  a  In  Eq.  (III-19)  may  be  used  to  determine  the 
velocity  of  sound  or  the  slowness  of  the  longitudinal  wave  In  the 
hexagonal  crystals.  Because  the  a  Is  also  the  same  for  different 
angle  the  velocity  of  sound  or  the  slowness  versus  the  angle  4> 
also  should  be  a  circle. 

Table  111-1  describes  the  relation  of  a  with  angle  e  at  some 
specific  direction  along  which  pure  longitudinal  waves  can  propagate 
In  hexagonal  crystals. 

As  soon  as  the  directions  of  the  pure  longitudinal  mode  are 
determined.  It  becomes  possible  to  derive  the  nonlinear  longitudinal 
wave  equation  whose  solution  describes  growth  of  the  second  harmonic 
of  a  pure  longitudinal  wave.  The  TOE  constants  can  be  obtained  from 
the  experiment  described  by  the  solution  to  the  nonlinear  wave 
equation. 

Suppose  we  have  ^  *  0  for  the  specific  directions.  Thus 
under  the  boundary  condition  of  the  pure  longitudinal  vibration  at 
the  surface  of  the  transducer,  that  Is  Uj'  *  AjSlnwt,  V|  *  WJ  *  0  when 
a"  *  0,  we  may  get  the  solutions  of  Eq.  (111-18)  as: 


Uj*  ■  A^s1n(kLa"  -  wt) 


U II  —  ull 

V1  "l 


(II 1-20) 


C.  r— 

Here  k^  ■  — ,  CL  "A2-  Is  the  velocity  of  longitudinal  waves. 

Putting  Eqs.  (III-7)-(IIl-15)  and  (11-34)  into  (11-41)  and 
using  (111-16)  with  ( I 11-20) »  we  can  get  the  second  approximation 
equation  of  the  pure  longitudinal  mode. 


An  estimate  of  the  magnitude  of  the  apex  angle  can  be  made  for 
specific  crystals.  If  one  Ignores  the  piezoelectric  terms,  one  can 
make  such  an  estimate  for  CdS,  one  of  the  most  Interesting  hexagonal 
crystals.  The  SOE  constants  of  CdS  are  given  by  Landold-Bomsteln 
(3)  as 


C11  * 

87 

GPa 

C12  " 

-54.6 

GPa 

C13* 

47.5 

GPa 

C33  * 

94.1 

GPa 

C44  = 

14.9 

GPa 

So  we  can  numerically  figure  out  from  (III-19),  when  e  *  ep  «  65.836°  or 
the  apex  angle  90°-65.836°  *  24.164°,  then  n&2  vanishes. 

ZnO  Is  another  Interesting  sample.  Its  SOE  constants  may  be 
found  In  the  same  reference  (3).  They  are: 


n  ■ 209 

GPa 

12  *  ,2° 

GPa 

IS’  104 

GPa 

33  ’  2,8 

GPa 

44  * 

GPa 

So  we  can  numerically  figure  out  from  (III-19),  when  e  «  ep  *  70.390° 
or  the  apex  angle  90°-70.390*  =  29.610°,  then  n$2  vanishes.  Figure  III-3 
and  Figure  II 1-4  show  the  coupling  coefficient  (e2/a)  versus  the  angle 
e  for  CdS  and  ZnO,  respectively,  in  each  figure  we  Indicate  the  pure 
mode  directions  along  which  02  vanishes. 


Figure  III-2.  Pure  longitudinal  mode  directions  in  hexagonal  crystals 
(indicated  by  arrows). 
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C 

C 

c 

c 

c 

c 

c 


111 

112 

113 

114 
124 
144 
222 


"  4C266  +  C112 

*  C122  “  Clll  +  C222 

*  2C366  +  C123 
“  2C156  "  3C124 

*  **C224  +  C1 14 

*  “2C456  +  C1 55 

*  4C166  +  2C111  +  C112 


( IV-2 ) 


If  we  substitute  (11-36)  Into  (IV-1),  we  would  be  able  to  give 
the  complete  expression  for  the  strain  energy  density  for  the  crystal 
of  trigonal  symmetry  In  the  twice-rotated,  or  double  primed,  coordinate 
system  as  was  given  In  the  previous  section  for  the  crystal  of 
hexagonal  symmetry.  But  we  already  know  that  if  our  only  concern  Is 
the  generation  of  the  second  harmonic  of  a  pure  longitudinal  mode,  it 
Is  sufficient  to  write  only  the  terms  containing  In  the  expression 
*3.  This  simplifies  the  expressions  a  great  deal. 

First,  let  us  look  at  the  expression  *2.  By  substituting 
(11-36)  Into  (IV-1),  the  expression  *2  can  be  given  In  the  double 
primed  coordinate  system  as: 

*  \  Cjj[ny2cos4e  +  (n^  +  ojjjjcos^ -2n^(n^'3  +  n3j)s1necos3e 
♦  (rjJ3  +  Ji3j)2cos2es1n2e]  +  j  Cj2(nJ2  +  n2*)cos2e 
+  CjjCn^cos^sIn^  +  +  n^coseslnefcos^  -  s1n2e) 

-  (n]|3  ♦  )2cos2es1n2e]  +  C^[n^'2s1n(2e)cos2es1n^(3cos24»  -  sin2*) 
+  +  n31)cos2e(cos2e  -  3sin2e)s1n*(3cos2*  -  sin2*) 
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2  7  7 

♦  3nJi(nJ2  +  n21)cos  es1necos4>(cos  <t>  -  3sin%) 

♦  \  (n^3  +  n3-|  )2cos(2e)sin{2e)s1n$(s1n24«  -  2cos2*) 

+  (nj2  +  n21)2cosesines1n<fi(s1n24>  -  3cos2$) 

+  (n^'2  +  n2l )(nl3  +  n3-j)cose(cos2e  -  2s1n2e)(cos2<fr  -  3s1n24>)cos4>] 

+  3-  C33[nJ2s1n4e  +  2nJ1(»jj3  +  n3^ )s1n3ecose 

+  (n^3  +  n31)2s1n2ecos2e]  +  C44[2n^,2cos2es1n2e 

+  (nj2  +  n22)s1n2e  +  2nJ1(nJ3  +  n31)cosesine 

x  (cos2e  -  s1n2e)  +  (nj3  +  n31)2(cos*e  +  s1n4e) 

-  4nJ3n51cos20s1n20]  .  (IV-3) 

The  strain  derivatives  can  be  calculated  and  simplified  by 
taking  into  account  n}'2  *  n21  and  n}'3  B  n31  as: 


A  O  p  9  9  9 

*  [C^cos  e  +  2C13cos  esin  e  +  2C14s1n(2e)cos  es1n«>(3cos%  -  slnS) 

♦  C33s1n*e  +  C44s1n2(2e)]nJ1  +  [-C11coses1n(2e) 

♦  Cj3cos(2e)s1n(2e)  +  2C^4cos2e(cos2e  -  3s1n2e)s1n<|> 

x  (3cos2$  -  sin2$)  +  C33s1n2esin(2e) +C44sin(4e)]n^3 

♦  [6cos2es1necos$(cos2$  -  3s1n2*)]nJ2  t  (IV-4) 


[3C14cos  es1necos<j»(cos  ♦  -  3s1n  +  [C^cos  e  -  C^cos  e 


+  4C^coses1nes1n<}»(s1nt^  -  3cost’$)  +  2C44s1n2e]n^2 

7  2  2  2 

+  [2C,,cose(cos  e  -  2s1n  e)cos4>(cos  ♦  -  3sin  ♦)]nV,  ,  ( IV-5) 

C-C11cos3esine  +  \  C13s1n(2e)cos(2e)  +  C^cos2e(cos2e  -  3s1n2e) 

•  s1n«t>(3cos2$  -  s1n24>)  +  C33s1n3ecose  +  C^sln (2e )cos (2e )3n^’-| 

2  2  2  2 

+  [2Cj4cose(cos  e  -  2s1n  e)cos$(cos  ♦  -  3s1n  ♦)]n^2 

+  [2C11cos2esin2e-  4C13cos2es1n2e  +  2C14cos(2e)s1n(2e) 

•  (sin2*  -  2cosZ*)s1n*  ♦  2C33cos2es1n2e+  4C44 

•  (cos4e  +  s1n4e  -  s1n2ecos2e)]nj3  .  ( IV-6) 


Carrying  out  a  similar  procedure  as  Mas  done  for  hexagonal 
crystals,  we  get  a  set  of  first  approximation  wave  equations  for  propa 
gatlon  along  the  a"-ax1s  In  trigonal  crystals: 


32u;  32v;  a2«; 
PnUV  *  a  j  ^  o  +  Bp  j 
0  1  3aB*  c  3a"* 


PnV?  »  Vi 


'l  »V  32w; 

4+3,  - 5“+3, - y 


(IV-7) 


O  «  C11cos4e  +  ^  C13s1n2(2e)  ♦  2C14s1n(2e)cos2es1n*(3cos2*  -  sin2*) 

+  C33s1n4e  +  C44sin2(2e) 

Y1  *  2"  C^cos2e  -  C^2Cos2e  +  2C,4coses1nes1n*(sin2*  -  3cos2*) 

+  C44s1n2e 

y2  *  ^  C^sin2(2e)  -  ^  3s1n2 (2e )  +  C^4cos(2e)s1n(2e)(sin2*  -  2cos2*) 

.  sin*  +  1  C33s1n2(2e)  +  2C44(1  -  |  s1n2(2e)) 

6^  *  3C14cos2es1necos*(cos2*  -  3s1n2*)  (IV-8) 

e2  *  "  7  cncos2es1n(2e)  +  \  C13cos(2e)s1n(2e) 

2  2  2  2  2 
♦  C^^cosce(cosce  -  3sin  e)sfn*(3cos  *  -  sin  *) 

+  j  C33s1n2es1n(2e)  +  ^  C44s1n(4e) 

2  2  2  2 
e3  *  C14cose(cos  e  -  2s1n  e)cos*(cos  *  -  3s1n  *) 

From  the  above  expressions  we  can  see  that  there  are  two 
coupling  terms  (e-j  and  Bg)  existing  In  the  longitudinal  wave  equation 
[the  first  of  Eq.  (1V-8)].  This  means  that  longitudinal  waves  excited 
by  the  transducer  may  couple  to  the  transverse  waves  of  both  perpendi¬ 
cular  and  horizontal  polarization  rather  than  only  one  polarization 
as  It  does  In  the  hexagonal  crystals.  So  In  this  case  we  seek  directions 
along  which  both  of  the  coupling  terms  vanish  Instead  of  only  one  term 
as  In  the  last  section. 


i  ',i  '.'J  ■.'.'  ■■>'■  'i.  ■ » 


II  ■  l  ■  1  •  !  •  . 
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First,  let  us  look  at  what  happens  In  the  basal  plane.  When 
e  *  0,  we  have 

3^  *  0 

and  (IV-9) 

&2  *  C^^sl n4>(3cos^4>  -  sin^$)  . 

Further, 

&2  *  0,  when  ♦  *  n  y,  n  *  0,  1,  2,  ...  (IV-10) 

This  means  that  both  $1  and  e2  vanish  and  the  pure  longitudinal  mode 
can  propagate  in  trigonal  crystals  only  along  two  special 
directions,  ♦  *  0  and  »/3,  In  the  basal  plane,  not  all  directions 
like  In  hexagonal  crystals. 

Second,  we  can  find  another  direction  of  the  pure  longitudinal 
mode,  which  Is 

e  «  (2n  ♦  1)  £  n  «  0,  1,  2,  3  ... 

Obviously,  all  of  the  directions  of  the  pure  longitudinal  mode  we 
found  above  are  common  for  any  trigonal  crystals  without  exception 
no  matter  what  their  SOE  constants  are.  It  was  pointed  out  (6)  that, 
besides  the  above  directions.  In  the  b-c  plane,  i.e. ,  4>  >  */2,  one 
unique  direction  may  also  be  discovered  with  vanishing  ^  and  @2. 
Fortunately,  we  have  already  had  the  possibility  to  find  it 
quantitatively.  It  Is  clear  that  when  +  ■  */2,  B-j  vanishes  and  e2 
is  reduced  to 


s2  *  "  7  cjicos2esin(2e)  +  7  C13cos(2e)sin(2e) 

-  C14cos2e(cos2e  -  3s1n2e)  +  j  C33s1n2esin(2e) 

+  ^C44s1n(40)  .  (IV-11) 

Thus,  the  direction  we  seek  Is  dependent  on  the  magnitude  of  the  second 
order  elastic  constants  of  the  sample  being  measured. 

The  quartz  Is  a  sample  whose  nonlinear  mechanical  properties 
are  of  wide  Interest.  The  SOE  constants  can  be  found  (3)  for  quartz 

as: 


Cn  «  86.6  GPa 
C-J2  *  6.7  GPa 

C13  ■  12.6  GPa 
Cj4  *  -17.8  GPa 
C33  «  106.1  GPa 
C44  *  57.8  GPa  . 


(IV-12) 


The  direction  with  vanishing  Bg  can  **  determined  numerically 
In  the  b-c  plane  from  (IV-11)  as 


ep  -  20.444°  . 

Figure  IV-2  shows  the  diagram  of  (B2/a)  versus  engle  e  In  the 
b-c  plane  for  quartz.  Lithium  nlobate  might  be  another  Interesting 
sample.  Its  SOE  constants  can  be  written  (3)  as  follows: 
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Cn  -  202.0  GPa 
Cj2  *  55.0  GPa 
C]3  -  71.0  GPa 
C14  ■  8.3  GPa 


(IV-13) 


C33  «  242.0  GP 
C44  *  60.1  GPa 


The  angle  ep  with  vanishing  $ ^  Q^ves 

ep  -  69.812°  . 

Figure  IV-3  shows  the  diagram  of  (^A*)  versus  angle  e  In  the 
b-c  plane  for  lithium  nlobate.  Finally,  we  give  a  summary  In  Figure  IV-4 
which  shows  the  direction  of  pure  modes  In  trigonal  crystals  in  general. 

So  far  we  have  determined  all  of  the  directions  of  pure 
longitudinal  inodes  In  the  crystals  of  quartz  and  lithium  nlobate. 

In  the  direction  with  vanishing  8-|  and  the  solution  of  the  first 
approximation  wave  equation  can  be  worked  out  with  the  boundary  con¬ 
dition  when  a"  *0.  Uj  «  AjSln(wt)  can  be  written  In  the  same  form 
as  (III-20): 

Uj  «  A^s1n(«t  -  k^a") 

(IV-14) 

V,  •  W,  »  0  . 

Using  a  similar  procedure  as  In  the  last  section,  we  may 
get  the  second  approximation  wave  equation  of  the  pure  longitudinal 
mode  for  the  trigonal  crystals.  But,  as  we  already  have  seen 
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Figure  IV-3.  Coupling  coefficient  |e9/a| versus  angle  e  for  LINbCL 
(at  i  *  (2n+l)n/2).  c 
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C 


Figure  IV-4.  Pure  longitudinal  mode  directions  In  trigonal  crystals 
(indicated  by  arrow). 
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above,  because  of  vanishing  and  Wj,  only  the  terms  Involved 

In  n must  be  kept  in  the  expression  of  the  part  of  the  strain 

energy  density  in  the  double  primed  coordinate  system,  as  well  as 

only  the  strain  derivative  of  elastic  energy  -Hr-  must  be  used  in 

3r»n 

the  developing  of  the  second  approximation  wave  equation. 

Therefore,  we  need  to  write  the  expression  for  4*3  only  with 

3 

those  terms  involved  in  That  is: 

♦3  *  C111cos6ecos24(cos24  -  3sin2$)2riJ3 

i  4  2  3  5  3 

+  £  C113cos  esin  en^  -  C^cos  esinen^ 

♦  2C124s1necos5esin<|>(l  +  2cos{2<fr))nj3 

+  7  Ci33cos2os1n4©r}^  ♦  2C13^cos3esin3esin4(4cos24  *  ))nj[3 
-  C144cos4esin2ecos24sin24n]|3  +  C155cos4esin2e(2  -  ^  s1n(2$))n}'3 

♦  g  C222cos6esin2a(l  *  4cos24)2n][3 


♦  Z  C333s1n6enll  4  2C344cos2es1nSnJ3 
4  J  C^4cos3esin3es1n4(l  -  4cos24)nf3  . 

And,  its  corresponding  strain  derivative  is 

J  C^cosScos^cos2*  -  3s1n24)2nJ2 
3  4  2  2  S  2 

*  7  Cii3c°s  0  sin  enj^  -  3C114cos  esinenJ-j 

♦  f  C12^s1necos5esin4(l  +  2cos(2*))nJ2 


(IV-15) 


A  v  v!*A* v\- v  v  v.-  v  /'.-v  ■  • 


♦  j-  C^jcos^esln^en^  +  6C^cos^es1n^0s1n<>(4cos^^  -  l)nj'2 

4  2  2  2  2 

-  3C14^cos  esin  ecos S>s1n%n^ 

+  3C155cos4es1n2e(2  -  ^  sin(24>))nj'2  • 

♦  C222cos6esin24’(1  “  4cos24)2nj2 

+  j-  C333Sln6en^2  +  6C~44cos20s1n*0nj] 

♦  4C^44cos3es1n3es1n4i(l  -  4cos2$)nj2  .  (IV-16) 


(IV-19) 


Kj  ■  &  -3*2  *  Cj^cos®ecos2*(cos2$  -  3s1n2$)2 

4  2  5 

♦  3Cj13cos  esin  e  -  6C^4cos  eslne 

+  3C124s1necos5esin$(l  +  2cos(2$))  +  3C133cos2es1n4e 

♦  12C.j3ijCos3e$1n3es1n$(4cos2$  -  1)  -  6C.j44cos4es1n2ecos2$s1n2$ 

♦  6C155cos4es1n2e(2  -  J-s1n(2*))  +  C222cos6esfn2$(1  -  4cos2*)2 

+  C333s1n6e  +  l2C344cos2es1n4e  +  8C444cos3es1n3es1n$(l  -  4cos2e)  . 

(IV-20) 

The  solution  of  Eq.  (IV-18)  can  be  obtained  with  the  boundary 
condition  UJJ  *  0  where  a"  *  0,  as 

(k.  A,)2 

U2(a".t)  «  -(3K2  +  K3)  awcos2(h»t  -  kLa")  .  (IV-21) 

Thus,  the  solution  of  the  nonlinear  wave  equation  of  the  pure 
longitudinal  mode  can  be  expressed  as 


(IV-22) 


(IV-23) 


The  parameters  K2  and  *3  for  the  specific  directions  of  a  pure 
longitudinal  mode  In  trigonal  crystals  are  written  In  tabular  form  In 
Table  IV— 1 -  The  directions  specified  are  the  directions  along  which 


Table  IV-1.  The  and  K.j  parameters  for  trigonal  crystals  along  the  direction  of  the  pure 
longitudinal  mode. 
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harmonic  generation  experiments  can  be  expected  to  yield  useful 
Information  about  third-order  elastic  constants  in  trigonal  crystals. 

V.  CONCLUSION 

The  results  given  In  this  Technical  Report  can  be  used  to 
evaluate  the  third-order  elastic  constants  that  determine  the  ampli¬ 
tude  of  the  second  harmonic  wave  propagating  in  the  listed  pure  mode 
directions  in  crystals  of  hexagonal  and  trigonal  symmetries. 

In  addition,  we  may  point  out  that  the  theory  developed  in  the 
early  chapters  of  this  report  is  sufficiently  general  that  one  can 
use  it  to  study  the  nonlinear  behavior  of  crystals  of  any  symetry. 

One  even  can  consider  piezoelectric  crystals  by  using  the  appropriate 
form  of  the  strain  energy  $  that  describes  the  effect  of  piezoelectric 
coupling  as  well  as  nonlinear  effects.  As  one  goes  to  crystals  of 
lower  symmetry,  however,  he  must  be  prepared  for  the  equations  to 
become  more  complicated. 
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